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1 Introduction 

Renormalization Group functions, i.e. /3-functions and anomalous dimensions, are of great 
importance in quantum field theory. Recently most of these functions have been computed 
at three-loop accuracy in the Standard Model (SM). The results for the gauge couplings 
have first been derived in [1, 2] and have been confirmed independently in [3]. For the 
top- Yukawa coupling and the Higgs self-interaction A the QCD, top- Yukawa and Higgs 
contributions have been derived in [4]. The result for the top- Yukawa /3-function has been 
extended to include the electroweak and all third generation Yukawa couplings in [5] where 
also the /3-functions for bottom and r- Yukawa have been presented. The one-loop and 
two-loop results for all SM couplings have been known for a long time [6-20] as have been 
partial three-loop results [21-26]. Four-loop /3-functions are available for QCD [27, 28] and 
the purely scalar part of the SM [29-31]. In this paper we present the extension of our result 
for the Higgs self-interaction and the anomalous dimension of the Higgs mass parameter to 
include the electroweak and all third generation Yukawa couplings as well. Especially the 
/3-function for the Higgs self-coupling is interesting because of its close connection to the 
question of vacuum stability in the Standard Model. It has been shown that the stability of 
the SM vacuum up to some energy scale A is approximately equivalent to the requirement 
that the running coupling A(//) > for \i < A [32-34]. Many analyses of this question 
have been performed [4, 35-41] during the last years. The main uncertainty stems from 
the experimental error on the top mass followed by the uncertainty in a a . But future linear 
colliders could greatly reduce these uncertainties and a possible Higgs mass of about 126 
GeV [42, 43] allows for both scenarios, a stable and an unstable (or to be more precise 
metastable) SM vacuum, within the present experimental and theoretical errors. For this 
reason we think that the present work will reduce the theoretical uncertainty connected to 
the running of A even further than our previous calculation [4]. 

In the following section the setup and the technical details of the calculation are discussed. 
After that we present our results and some numerics in order to determine the significance 
of the new terms. 
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2 Calculation 



The gauge group of the SM is an SU C (3) x SU(2) x U y (l) which is spontaneously broken 
to an SU C (3) x Uq(1) at the Fermi scale. Our calculation is performed in the unbroken 
phase of the SM which is justified by the fact that the UV behaviour and therefore the 
renormalization constants for fields and vertices do not depend on masses in the MS-scheme 
[44]. 

The Lagrangian of the SM can be decomposed into the following pieces: 

£■ = C-QCD ~\~ C'EW ~\~ £-Y ukawa + (2.1) 

The QCD and electroweak (EW) part are implemented in the usual way with the gauge 
fields (SU C (3), a = 1, ... ,8), (SU(2), a = 1, 2, 3) and (U y (l)). These appear in 
the covariant derivative 

D» = ^ - igiY f B^ - i^a a W afl - ig s T a A af * (2.2) 

with the Pauli matrices a a and the hypercharge Yf of the field / on which the covariant 
derivative acts. In the Yukawa part we neglect the first two generations and the mixing of 
generations. The W-fermion-vertices are taken to be diagonal in the generations as well, 
i.e. we set the CKM matrix to the unit matrix. Light quarks and leptons are present in 
the QCD and electroweak sector however. This leads to the Lagrangian 

\ _ J 1 J (2.3) 

- y T \t r &L l + L L m R j 

for the Yukawa sector. The complex scalar field <1> and the left-handed quarks and leptons 
are doublets under SU(2): 

-(£)• ^(:)- <-» 

The indices L and R indicate the left- and right-handed part of the fields as obtained by 
the projectors 

^ = ^(1-75) P R = \{1+^). (2.5) 
Finally, we have the Higgs sector with 

= (D^)\D^<S>) _ m 2 $t$ _ x (V$) 2 . (2.6) 

For every field and vertex a counterterm is introduced and the corresponding renormal- 
ization constant is calculated order by order in perturbation theory. The left-handed and 
right-handed parts of fermion fields, quark-gluon- vertices and fermion-B- vertices are renor- 
malized with different counterterms. The renormalization constant for the gauge, Yukawa 
and Higgs couplings can be obtained in different ways, e.g. 

z (ttB) z (ttB) 

91 ~ 7 (2r) fZ(^Bj ~ 7 (2t) JZ^B) ~ ' ' ' 
Z 2,L V ^3 Z 2,_RV Z '3 
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where zf 1 "' : ' n is the renormalization constant for the vertex of the (renormalized) fields 

(2 f) 

fx, . . . , f n , L j R the field strength renormalization constant for the left-handed (L) / right- 
handed (R) part of the fermion field / and Z^ 9 ^ the field strength renormalization constant 
for the gauge field g. 

Likewise, the renormalization constant for Yukawa couplings can be computed from the 
renormalization constant for any vertex proportional to this coupling and the renormaliza- 
tion constants for the external legs of this vertex, e.g. 

\l ^2,L A 2,R ^2 

where Z^^ is the field strength renormalization constant for the scalar doublet. All 
renormalization constants are defined in a minimal way as 

Z = 1 + 6Z, (2.9) 

with 5Z containing only poles in the regulating parameter e = (4 — D)/2 where D is 
the engineering space-time dimension. The Higgs self-coupling A is renormalized with a 
counterterm 5Z\ that is not proportional to A but also has terms proportional to four 
Yukawa couplings. Consequently, this is a feature of the corresponding /3-function as well. 
We find 

-X B K$ B ) 2 =(-A + «5Z 1 (4 * ) )($t^ 2 



=► -(X + SZ X ) (z™) =(-\ + 5z{^) (2.10) 
^^ = (f»)" 2 (A-< $) )-A, 

with the counterterm 5Z^^ for the four-<3?-vertices and the index B marking bare quan- 
tities. The /3-function for any coupling X is defined as 

and is given as a power series in all considered couplings of the SM, i.e. g s , gi, gi,y t ,yb,y T 
and A. The mass parameter m? of the scalar field is neglected in the calculation as it has 
no influence on the UV behaviour of the couplings in the MS-scheme. 

The /3-function (or anomalous dimension) describing the running of this mass parameter 
m 2 in eq. (2.6) can be computed from the renormalization constant of the local operator 
02$ := An insertion of 02$ into a Green's function, e.g. with two external ^-fields, 

is renormalized as [C>2$] = Z^02^ where [O2*] is the corresponding finite operator. From 
[0 2 «] = Z m 20^ c and 0^| rc = zf <s>) 2 <s> it follows that 



Z m 2 = (z { ™ ) ) X Z&. (2.12) 
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For this project we need the renormalization constants 5Z\ , Z^i and at three-loop 

level and all other counterterms at two-loop accuracy or less. We perform our calculation 
in a general i?£-gauge with different gauge parameters £1, £2 and £ for the gauge fields B, 
W and A. The /3-functions for all couplings are independent of the gauge parameters which 
serves as an important check for the result. In order to compute all counterterms up to 
two-loop level and the one and two-loop diagrams with counterterm insertions contributing 
to the three-loop result it has been convenient to use a setup where different isospin con- 
figurations of a field and different fermion generations (in our case the distinction between 
third generation and light is enough) are implemented as separate fields as many countert- 
erms depend on those. Additionally, the U y (l) hypercharge depends on the isospin and 
differs for the left-handed and right-handed part of the field. So we use the set of fields 

UM, e-,Ve,T,v T , A a ^B^W^,W 2fl ,W^ , c a , c^c^cj, , $ a , $ 2 

quarks lcptons gauge bosons ghosts 

and their anti- fields. The fermion fields have to be split up in a left-handed and a right- 
handed part during the calculation. The price we pay, however, is that many diagrams are 
produced which look the same in momentum space if no counterterms. are inserted. For 
the 1PI process with four external <&-legs ~ 2.3 x 10 6 dia grams are generated at three-loop 
level. In order to reduce this number and because we do not need three-loop diagrams with 
counterterm insertions anyway we have chosen a second smaller set of fields for this part 
of the calculation 

^ , , A a »,B»,W a » ,c a ,c? v ,<S> i . 

quarks leptons gauge bosons ghosts 

Here the index i marks the isospin of <I> and the left-handed fermions. For the right-handed 
fermions i = 1, 2 is just a label to mark the flavour, e.g. q\ t R = tn, (72, i? — &R- Fermion loops 
without Yukawa interactions are multiplied by the number of generations N g in order to 
include the light fermions. The indices of external particles can be explicitly chosen. Using 
this setup only 573692 diagrams are produced for the four-$i-process at three loops. The 
computation of the SU(2) x U y (l) group factors has been implemented with Mathematica 
using labels for the left-handed and right-handed part at each quark- Yukawa-vertex and 
.B-fermion-vertex as well as for the three different structures in the four-VF-vertex. With 
the help of these labels we can completely factorize the SU(2) x U y (l) part from the 
momentum space diagram. The QCD colour factors have been calculated with the FORM 
package COLOR [45]. All Feynman diagrams have been automatically generated with 
QGRAF [46]. 

In order to check our setup we have computed all Yukawa and gauge coupling renormaliza- 
tion constants at two-loop level from at least two different vertices with the first set of fields 
and compared the result to the literature. The same has been done for the renormalization 
constants of the gauge, ghost and scalar fields. We also explicitly checked that we get 
the same renormalization constants for $1 and $2 as well as for the left-handed fermion 
flavours of the same generation at two loops. Another check has been the finiteness of the 
three- £>-vertex up to two loops. 
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The renormalization constant for $ has been computed with both sets of fields at three- 
loop level which yields the same result as in [2]. The renormalization constants for the 
W, cw-i 3*1 and $2 fields as well as for the W-cvy-cvy-vertex, the C^-^i-^i-vertex, the 
C^-^i-^i-vertex and the 4-<3?i-vertex have been computed with both sets of fields up to 
two loops with the same result. 

As explained in detail in [4] some diagrams with four external ^-fields where two external 
momenta are set to zero suffer from IR divergences which mix with the UV ones in dimen- 
sional regularization. We therefore use the same method as in [4] and introduce the same 
auxiliary mass parameter M 2 in every propagator denominator. Subdivergences oc M 2 are 
canceled by counterterms 

^SZ^ —5zfP WZW^, —bZ^ B„B» and —bZ^ (2.13) 

9 m z m ' 9 M " 2 2 

Counterterms oc M that would arise for fermions cannot appear because there are no M in 
the numerators of propagators. Ghost mass terms ^f-fi^^a c a c a for the SU(3) and SU(2) 
ghosts do not appear because of the momentum dependence of the ghost-gauge boson- 
vertex. The remaining divergences are the mass-independent UV ones we are looking for. 
This method has been suggested in [47] and has been elaborated on in the context of three- 
loop calculations in [48]. The resulting massive tadpole integrals can be computed with 
the FORM-based program MATAD [49]. 
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Figure 1: Some diagrams contributing to the renormalization of the 4-<l>-vertices 

As opposed to the case of the Yukawa coupling /3-functions (see [4, 5]) a completely naive 
treatment of 75 in dimensional regularization is possible for f3\ and f3 m 2 . In four dimensions 
we define 

75 = i7°7 1 7 2 7 3 = jf»vpvYl v Yl a with e i23 = 1 = ~e m2 \ (2.14) 

In order to have a non-naive contribution from a fermion loop with a 75 matrix in it at 
least four free Lorentz indices or momenta on the external lines of the minimal subgraph 
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containing this fermion loop are required. These can be indices from the gauge boson 
vertices or the internal momenta from other loops which act as external momenta to the 
minimal subgraph containing the fermion loop in question. External momenta of the 
whole diagram can be set to zero as the renormalization constants in the MS-scheme do 
not depend on those. In four dimensions the trace of such a fermion line will produce a 
result oc £^1^2/^4 wnere Ml) • • • > A*4 are the aforementioned free Lorentz indices. If there is 
a second fermion line the trace of which also yields an e-tensor these two e-tensors can be 
contracted and we can get a non- naive contribution from 75. 

Let us consider two examples. The diagram in Fig.l (c) has two fermion loops. If we take 
one of the fermion loops with the momenta on the two external $-legs set to zero we have 
two indices from the gluon lines attached to the fermion loop and one loop momentum going 
through the two gluons and acting as an external momentum to the subgraph containing 
only this fermion loop. This is not enough to have a non- naive 75 contribution from this 
graph. 

The fermion loop in Fig.l (f) has three Lorentz indices and two external loop momenta 
but there is no second fermion line to produce a second e-tensor. Furthermore - as we 
set all external momenta to zero - there are no free Lorentz indices or momenta in the 
final result to support an e-tensor there. Therefore any contribution with an e-tensor from 
this fermion loop (which is the only antisymmetric Lorentz structure) must vanish after 
contraction with the Lorentz structures from the W-bosons. 



3 Results 

In this section we give the results for the three-loop /3-functions for the couplings A and 
the mass parameter m 2 setting all gauge group factors to their SM values. All results of 
this work can found at 

http:/ /www-ttp. particle. uni-karlsruhe.de/Progdata/ttpl3/ttpl3-008/. 1 

We denote the number of generations by N g . 



1 There we also present the results in a form where the QCD colour factors are not set to numbers but 
generically expressed through the quadratic Casimir operators C F and Ca of the quark and the adjoint 
representation of the corresponding Lie algebra, the dimension of the quark representation d R and the trace 
T F defined through T F 8 ab = TV (T a T b ) with the group generators T a of the quark representation. 
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A W = - vt - vt* + 9ty Q + 9l9l\ + gf^ + Ay 2 2 + Ay 2 6 - Xgf- - Xgf- 

+ A 2 12 + y t 2 A6-y 4 3, 
# =y^ + y ft 6 15 - givll - glyf- + gl(^-- 2N a ) - + gfyf- 



11 ooo9 n vt / 9 1 2 \ j i}25 /io5 



, 222 ii , 2 2 2 y 2 4 f y ' , Z at \ 42^°, 42° 

4 2 /239 10 \ 6 /59 10 \ x 4 1 , 4 3 
+ Ayfy 2 ^ + A 5 iy 2 ^ + \g% ( + 5iV s ) + A 5l 2 y 2 ^ + Ay 2 ^ 2 ' 



4 4 16 V 4 aiab 12 

+ Xgfgfj- + Ay 4 + |jV.) - A 2 y 2 24 - A 2 y 2 72 + A 2 y 2 54 

+ A 2 y 2 18 - A 3 156 - y t 2 y 4 3 - y*gf- + y 2 ^ - y^gff 

- y 2 Ay 2 21 + y 2 Ay 2 2 f + - y t 2 A 2 72 - y 4 y 2 3 - y\g\\ 

- vt>\ + yfl5 - y 2 y 4 16 + ff 2 Ay 2 40 + g*y?\40 - y 2 y t 4 16. 



*( 143 1rt> \ 2 fi 297 4 . fi 2 297 „/ 1599 „ A 
- l2Csj - ^Y 6 -^- - vtvt™ - y e b y 2 T — + y 8 b ( — — - 36C 3 

2fi /ll37 A \ 2 fi /34ll „ \ 4 4 /4503 273 ^ 13 Ar 

+ s 2 Y 6 (-32- " 9C 3 ) + g\yl - 27 C3 ) + giyt (— - — Cs - T iV s 

4 2 2 9 4 4 /13653 819 A 39 A7 \ B 2 ) 5739 99 9 Ar 

+ <? 2 YY 2 g + sk 4 (^g- - ^Cs - T iV 9 j + g\yl (- w + T Cs + ^ 

62 / 17217 297 27 \ B / 982291 2781 14749 AT 

+ 9 ^ {~m + T Cs + +52 1^072- " 128 <3 - ~m~ N ' 

-45A^3 - § JV?) + g\£ + 33C 3 ) + gU - 25 C 3 

2 2 4 / 1" 381 \ 2 2 2 2^ 224/ 3239 311 

+ 5i5 2 y T l~^~~^ 3 ) - 9ig2V b y T ^ + 9i92V b 

2 4 2 ^1833 3 A 1 A \ 2 4 2 ^4179 5 Ar 
+ Sisll/? h^T - oCs - ^N g + y 2 9 4 y 2 — r + 9C 3 + -N g 



(3.1) 



V 256 2^° 2 V 3i ^ 6 V 256 "° 2 

26 / 54053 405 8341 10 Ar2 \ 4 4 Z5697 375 , 65 

+ 9 ^ {- M56" " ^ Cs " "864 ^ "27^0 + ^ + TtT C3 + 12 

4 2 2 4 1 4 4 /15137 2035 ^ 415 Ar \ 4 2 2 Z6657 15 . 5 Ar N 
+ ("3456 - " W Cs - ^iV s J + 5l <7 2 y T ( W " yCa - giV^ 

/4403 9 A 25 \ 44/ 64693 873 , 149 Ar ^ AT A 

50 ,.2\ 6 2/ 3y 29 15 55 A \ 6 2 / 12043 5 A 95 A7 

tA { 29779 75 18001 61 MO^, 
6912 32 Sd 2592 9 9 sS 243 9 



o / 6845 99 A 20735 95 . 125 Ar9 \ x fi / 1241 

+ 9i + C 3 N. + — iV C 3 N 2 + Ay 6 66C3 

yi V 9216 128 s 1728 9 9 9S 81 V y " V 8 s 

+ A^Y 240 + Ay 4 y T 2 240 + Ay 6 6 - 198C 3 ) + Ay 2 V 4 (-^ + 171(a) 

2 2 2 v , x 2 4 / 4977 , \ , x 4 2 ( 1311 117 . 21 



- Xgiyivt 27 + A^y? — + 513C 3 ) + (--^ -Cs - -jN, 

x 4 o / 3933 351 63 \ . fi / 46489 2259 A 3515 AT 

+90iY sC 3 + yiV 9 2 ) + Ay 2 y 4 - 1170,) - A 5l 2 y 2 y 2 9 + Ay 2 y 4 (-^ 
222/ 3771 , iocA \ , ,222/ 3009 



+249C 3 ) + Xghhr ( ~ + 126C3J + \g\g\yl f — ^- + 12C3 

, 2 4 /4553 249 ^ 33 AT \ ,42/ 1783 123^ 6 5 "\ 

yii/6 V 1728 6 s 36 V yiy H 8 8 s 4 9 9 ^ 
, fi /12679 9 A 5995 AT 190 , T , 1750 n A l2 4 /717 \ 

- A 2 y 2 y T 2 216 + A 2 y 4 + 756 C3 ) + " 144C 3 ) + A 2 9 2 y 2 

-432C 3 ) + A 2 y 4 - 513Cs - 141JV.) + A 2 9 2 y 2 (-^ + 96C 3 

t 2„2„,2 / 417 1noA \ , \2„2„2 l qqq ictaI i \2„4 



+ A 2 <? 2 y 2 ^— - 192C 3 J + A 2 y 2 <? 2 (-333 - 162C 3 ) + A 2 y 4 (-183 - 81Cs 
235 \ 

~~3~ N9 J + A3y ' 291 + A3yt2873 + 92 ( " 474 + 72Ca) + ^ ( " 158 + 24Cs) 
+ A 4 (3588 + 2016C 3 ) + M (-^) + y t 2 yfy 4 12 + y t 2 y 4 y T 2 f + y t 2 y 6 6 (- ^ 

2 2 4477 2 4 2V 2 4 2 ^ 351 117 



-36C 3 ) + yigivt— + yfgh;^ + vtaht + ~ Ca - 

2H f 17217 2 9 7 27 \ 224/ 2299 \ 

+ ^ 2 (-^56- + — Cs + T 7 + ^ (~9T + 26C3 ) 

, 2 2 2 2 29 2 2 2 2 /1001 31 \ 224 /3103 27 1 \ 
+ vt9i92V;-£ + V t 9i92Vt y-QQ- + yCs J +y t 9i92 [^r + ^ 3 + 2 9 J 

4_ 2 4 2 701 2 4 2 f 709 j. A 2 4 2 /23521 5 

+ ^ 2ff4 ( Hf - ^ 3 + ^1^) + y * Xy >° + y2M21 

+ y 2 M i^f- + 144C 3 ) - y 2 \g 2 2 y 2 27 + (-^ + 54C 3 

2, 4/ 3933 351 , 63 AT \ 2x 2 2 « 2 , 2 2 / 929 „. 
+ 2/ t 2 Aff 2 4 ^Cs - ~rN g - y 2 \g 2 y 2 T 9 + y 2 A 5 2 y 2 -— - 2C 3 



64 2 * 4 9 7 ' Jt V 12 

2x2 2/ 6509 2l 4 / 112447 449 635 Ar 

+ « (-^- + 177C 3 ) + y f 2 Ay 4 - -Cs - -^N g 

- y 2 X 2 y 2 216 + y t 2 A 2 y 2 (117 - 864C 3 ) + y 2 A 2 y 2 2 - 432C; 
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+ V 2 A 2 9l [~ - 48Cs) + y 2 A 3 873 - ^72 + yfy 2 b y 2 f + yfy^ 

4 2 2 477 4 4 /13653 819 , 39 \ 4 , , /1337 \ 
+ vUlvl^ + Vhi - -^Cs - -jN g J + yfgfy 2 - 28C 3 J 

4 2 2/ 1079 743 \ 4 4 /100913 2957 A 115 „ 
+ V t 9i92 ~ + V ^ [-M56- + W C3 " 36 

+ ^ r 2 240 + yf Xyl {^f- + 144(a) + v'M [~- + 513 C 3 
+ vW, (-^f + 57 C3 ) + , t 4 A 2 (™ + 756 C3 ) - M*Z 
+ « ("f " 36Ca) + rfrf - 27 C3 ) + rfrf + 17 6 ) 

+ W?A f ^ - 198C 3 ) + y t 8 f- ^ - 36C 3 ") + rfrf (-38 + 240C 3 ) 



+ sf.sivt {—^ + — Cs ) + gtnttivt [— - 36Cs ) + 9:9(92 {-y n ° + 6iY * C3 

2 4 o /683 \ , 4 3 / 187 Ar 22 AT \ 2 fi / 187 „ 

+ 9 2 Avl [-24 " ISCsJ + ff s 2 5i5l (-^4-^ + y^CaJ + <fo! (-^4^. 

22 \ / 489 \ 

+yA r 9 C 3 ) + fl?Ay? (895 - 1296C 3 ) + g 2 s \g\vl [-— + 216( 3 ) 

^ 4 /135 _ \ 2 A „2, /2 A_991 \ 2 4 f5j> 



+ rfflirf ( -| + 24C 3 ) + ^ - 54 C3 ) + (fcg (-^JV. + 18iY 9 C 3 

641 | 136 ^ ^ | ^2„2„2„.2 /233 Qft/ ._"\ , „2„2„4 ' "1 
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+ £A<£ iV 9 - 72iV 9 C 3 J + ^Atftf + 40C 3 J + g z s \gt [jN g 

-f iV 9 Cs) + 5 2 A 2 y 2 (-1224 + 1152Cs) + (-2 - 48C 3 ) 

+ glvhlvl (-8 + 96C 3 ) + gtvhi {^f- - 54C 3 ) + <7 2 y 2 <? 2 <? 2 2 - 36C3 
+ glvh\ - I8C3) + g 2 y 2 M (82 - 96C 3 ) + g 2 y 2 \g 2 2 (-^ + 216C3 

+ g 2 y 2 Xg 2 ("^ + 136Cs) + <? S \ 2 A 2 (-1224 + 1152Cs) + glvtvl (-2 - 48Cs 

+ alvUl (-§ + 24 Cs ) + g'yfgl - f Cs) + rfrfA (895 - 1296C 3 ) 

+ ^ (-38 + 240C 3 ) + gtvt [~ + 32C 3 + 40iV 9 ) + g^Xy 2 (*|° 

/ 1 82(1 

-48C 3 - 64iV 9 ) + 0^192 + 5s 4 y 2 A ^— - 48C 3 - 64iV 5 
+ 5 4 y 4 (-^f + 32C 3 + 40iV 9 ). (3.2) 

The purely A-dependent part of eq. (3.2) has been derived before in [29, 30], the full 
one-loop and two-loop results in eq. (3.1) are in agreement with [18-20] and in the limit 
<7i) 52) Vb, y T — > we reproduce the result [4] at three- loop level. 
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The running of the m 2 parameter is given by 
8 {l) Q 1 

H=yl + y^-9l- A -gl- A + ^ + y% 

m z 4 4 

P ( 2 4 9 4 27, 2 2 15 2 2 45 

a ( 385 5 \ 2 2 2 2 2 2 

+ 52 (,-"32" + 2^J + ™r ¥ + 24 + 5152 16 (3.3) 

/ IK7 25 \ 

+ Ql V^6 + T8 Nb ) ~ Ay ' 12 " Xy2bM + A5236 + Xg2lU 
~ A 2 30 - + y y 2 f + rfrfg - y 2 A36 
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(3) 

^ = y " (-if + i5Ca ) + y2 ^ 72 + M72 + J + 45Cs ) 

+ glvt \~ + 54C 3 ) - <?kY 2 f + alvt [~ + 162C3 

4 2 / 255 81 21 \ 4 2 / 765 243 , 63 



+ 92V T [~ m ~ T Cs - Y N g j + g 2 y b - — Cs - ^N g 

a( 39415 711 2867 Ar , rArA 35 Ar2 \ 2 4 /291 

2 2 2 9 2 4 / 1067 _ \ 2 2 2/^ 2331 , r ,\ 2 2 2 865 

- sfy 2 ^ + 9lVb ("IF + C V + 9i92Vt + C V " 9l92Vb ~u 

2 4 /2691 405 ^ 21 , r „\ * of 3607 15 > 6 5 

4 2 / 79207 35 . 155 , \ 4 2 /1053 207^ 55 , r 
+ ^ {- " l2^ 3 " "72-^0 + ^ 2 ("ST " T^ 3 + ~8 9 ~ 3N & 

c / 839 51 4375 AT 95 , T A 875 »\ , 4 /261 „ n , 
+ 9 "{W8- 16 Cs + " Y^ C3 + 243*?) + ^ I" + ?2CS 

r^i + 4- ^ m - 4 1^ r 567 



- A^^108 + Ay* ^— + 216C 3 J + \givi ^— - IO8C3 J + Xgivt ~ 32 4C 3 

, 4 / 11511 A 153 \ , 2 2 / 549 \ 4 2 2 /393 

+ X9l [~^- ~ 162Cs - — N g j + Xgjy 2 (-- - + 36CaJ + Ay 2 y 2 (^— - 132£ 

, 2 2 / 1701 \ x 4 / 1077 A 85, \ , 2 2 99 l2 2 297 

+ A 5 ?pl (-^g- + 36C 3 J + Ay? (-- — - I8C3 - yiV 9 J + X 2 y 2 r - + X 2 y 2 — 

+ X 2 g 2 2 (-63 - IO8C3) + X 2 g 2 (-21 - 36C 3 ) + A 3 1026 + y 2 y r 4 72 + y t 2 y 2 y 2 y 
2 4 f 4047 \ 2 2 2 / 27\ 2 2 2 / 243 

2 4/ 765 243 A 63\ 2 2 2 9 2 2 2 /" 929 > 
+ 2/*52 I - - — C3 - -g-iV 9 1 - y t giy T - + y t 9iy b I -y- - C3 

222/ 3277 117 \ 24/ 123103 149 A 635 , \ 2 , 2 
+ . 2 5l 2 9 l (-^ + — CaJ + y t 2 ^ (- ^ - ^C3 - ^ 9 j - y 2 Ay 2 108 

+ y 2 Ay 2 (-^ - 2I6C3) + y t 2 Ay 2 - 324 C3 ) + y t 2 Ay 2 (-^ - 6OC3 
.2x2 297 , , „ 4. 2/4047 , Q _\ , , 2 / 3177 



y'A 2 — + y?y T 2 72 + y y b + 36C 3 J + ytd [~^- + 162C 3 J 



+ M {-^ + 12Cs) + itf A + 216Ca) + V? + 45 C 3 

+ aht - 360C 3 ) + glg\yl (-^ + IO8C3) + y 2 y 2 4 (^iV a - 36iY 9 C 3 ) 
+ iahl i~ + 2 0C 3 ) + 9 2 A (f A. - f iV 9 C 3 ) + (-612 + 576C 3 ) 

+ ^ (41 - 48C 3 ) + y 2 y t 2 y 2 2 (-^ + IO8C3) + ghhl [~ + 68C3 

+ y 2 y t 2 A (-612 + 576C 3 ) + gW t - 360C 3 ) 



I (*f - 24C 3 - 32iV 9 ) + 9ty 2 (j£ ~ 24Cs " 32iY « 



(3.4) 
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The one- loop and two-loop parts of this result are in agreement with [19, 20]. The purely 
A-dependent part can be found in [29, 30]. and for g\, g2, y b , y T — > we reproduce the result 
[4] again. 

Now we want to give a numerical evaluation of the /3-functions at the scale of the top mass 
in order to get an idea of the size of the new terms. For M t ~ 173.5 GeV, Mh ~ 126 GeV 
and a s = 0.1184 [50] we get the couplings in the MS-scheme at this scale using one- loop 
matching relations [51-53]: 2 y t (M t ) w 0.94, g s {M t ) » 1.16, g 2 (M t ) « 0.65, gi(M t ) « 0.36 
and A(Mf) ~ 0.13. The lighter Yukawa couplings can be estimated from the MS-masses 
m b » 4.18 GeV and m T « 1.777 GeV [50] to be y b « « 0.02 and y T « ^/2^ « 0.01 

For = Mt) we find at one-loop order 



^\lM = M t ) 

(16tt 2 ) 



-1.5 x 10" 


" 2 +4.4 x 10 


_3 -1.6 x 10' 


~ 3 +1.3 x 10~ 3 


% 4 




9§A 


A 2 


+6.3 x 10" 


" 4 -1.6 x 10 


" 4 +1.3 x 10' 


~ 4 2. x 10~ 5 




A3? 


2 2 

9291 


9 4 


+2. x 10~ ( 


! +1.6 x 10" 


7 -3. x 10~ 9 


-6.3 x 10~ n . 




Ay? 


^ 4 





(3.5) 



At two-loop order the largest terms with y b are of O(10 7 ) and those with y T of O(10 9 ). 
Neglecting these small terms we find 



(16vr 2 ) 2 



-6.8 x 10 


_4 +4.1 x 10" 


" 4 +2.5 x 10" 


4 -4.3 x 10" 


-5 


1 4 




9lH 


A 2 2/ 2 




+2.2 x 10" 


~ 5 +1.5 x 10" 


" 5 -1.4 x 10" 


5 +1.1 x 10" 


-5 


9^9 1 


9 2 2 A 2 




9% 




+1. x 10" 


5 +9.6 x 10"' 


3 -7. x 10" 6 - 


-6.1 x 10~ 6 






9|A 


v ^ 






-5.4 x 10 


" 6 -4.6 x 10" 


" 6 +4.2 x 10" 


6 -2.6 x 10" 


-6 


4 2 


9 4 9? 


Hoi 


9\9\ 




+ 1.9 x 10 


_6 +1.6 x 10" 


" 6 -1.4 x 10" 


6 +1.4 x 10" 


-6 




A 2 9 2 


vh\ 


9|Ag 2 




-6.2 x 10 

" v 


7 + (terms oc y 2 ,y 2 ) . 







(3.6) 



9? 

At three-loop level we only give the largest terms and omit small ones of O(10~ 7 ): 



2 This has been done with the Mathematica package SMPoleMatching.m by F. Bezrukov which can be 
downloaded at http://www.inr.ac.ru/~fedor/SM/download.php 



- 12 - 



(3) 

^ (16tt~2)3^ = 5.9 x 10~ 5 ^-3.8 x 10"^ +2.5 x 10^-2.3 X 10~ 5 , 
givt vl aivi g^vi 

+8.8 x 1(T 6 +5.9 x 10~ 6 # +5.5 x 10~ 6 -5.1 x 1(T 6 

-4.4 x 10-^+4.2 x 10~ 6 -1.9 x 10~ 6 +1.7 x 10~ 6 (3 7) 

9§ 9|A g\\y\ g\y\ 

+1.5 x 10~ 6 # +1.3 x 10~ 6 -1.3 x 10~ 6 +1. x 1(T 6 

9lkvl ~ST alvUl gtvtst 
-9. x 1CT 7 -9.2 x 10~~ 7 +8.9 x 10~ 7 +8.3 x 10^ +smaller terms . 

9|9?A g* A 2 929?!/? 9? A2 ?/ 2 

The dominant contributions contain only g a , y t and A as suggested in [4] . This result there- 
fore explicitly confirms the validity of the approximation made in [4] at the scale of the top 
mass. On the other hand some terms proportional to powers of g\ are of the same numerical 
significance as some of the electroweak two-loop contributions and in including them we re- 
duce the theoretical uncertainty in the evolution of A stemming from f5\ even further. The 

overall value of the three-loop /3-function at this scale ^p^gjr^ = 4.2 x 1(T 5 is about 17% 

= 3.6 x 10~ 5 ). 

92,91^0 

The whole /3-function is dominated by the negative one-loop term oc yf but with the new 
terms presented here A will decrease less which means that these contributions slightly 
enhance the stability of the electroweak vacuum state in the SM. 

For /3 m 2(fi = Mf) we find at one-loop order 

/3 (1 2 V = Mt) „ „ „ 

1.7 x 10~ 2 -6. x 10" 3 +4.9 x 10~ 3 -6.1 x 10~ 4 



(3) 

larger than with the electroweak contributions neglected ( ^n^^T^ 



(167r 2 )m 2 

Vt ai x 9? (3.8) 



+7.6 x l(T 6 +6.3 x 10 -7 



Omitting the small contributions proportional to y b and y T the two-loop terms look like 
this: 

P (2) 2 (/i = M t ) 
/ on 9 9 = 9-6 x 1(T 4 -2.1 x 10" 4 -1.7 x 1(T 4 
(16-K 2 ) 2 m 2 v v v v ' 

9 2 y 2 vt ^v 2 

+8.4 x 10" 5 +7.9 x 10~ 5 +2.1 x 10~ 5 +1.6 x 10~ 5 -2. x 10~ f 

9 2 Vt 9 2 2 a at vhl * 2 

+8. x 10" 6 +6.6 x 10" 6 +2. x 1(T 6 +(terms oc y 2 b ,y 2 T ) . 
Ag 2 gf g%g 2 



(3.9) 
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Again we only give the largest terms at three-loop level and omit small ones of O(10 7 ): 

3 i3 )(u = M+) , - . . 

= -5.6 x 10" 5 +3.4 x 10~ 5 +2.7 x 10~ 5 +1.3 x 1(T 5 



(167r 2 ) 3 m 2 



9 2 s yf 9iv 2 t v b t gl 



+9. x 10~ 6 ,-5. x 1(T 6 -3.9 x 10~ 6 -3.5 x 1(T 6 +3.2 x 10~ 6 -1.7 x 10" 

aial #M ~aP~ 

+9.7 x 10~ 7 +5.7 x 10^ 7 +5.7 x 10~ 7 +5.6 x 10~ 7 +smaller terms . 



9%9 2 Vt A3 alvhl 



(3.10) 



We see that especially the term oc g\ is not much smaller than the three dominant terms 
that have already been computed in [4]. The overall three-loop results with gi-,9\ ~ > 
and the electroweak interaction switched on differ only by 1. x 10~ 6 due to cancellations 
between the new terms. 



4 Conclusions 



We have computed the three-loop /3-functions for the quartic Higgs self-coupling and for the 
mass parameter m? in the unbroken phase of the SM, neglecting only the Yukawa couplings 
of the first two generations and the mixing of quark generations (an extension that could 
easily be made but which is numerically negligible). The electroweak contributions are 
small as expected which confirms the validity of the approximation made in our previous 
calculation [4]. The full three-loop /3-function for the quartic Higgs self- interaction is about 
17% larger than without the electroweak terms thus making the SM vacuum state slightly 
more stable. 

Note added: Recently, a similar calculation has independently confirmed the results pre- 
sented in this paper [54]. 
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